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1. INTRODUCTION 

In recent years, research on self-balancing two-wheeled bicycle has been interested by many 
scientists. In particular, a difficult problem is the study of self-balancing problem of the robot. To solve the 
problem of balancing two-wheeled bicycle, there are three basic methods as follows; (a) controlling balance 
by the flywheel, as in the studies of Beznos [1], Xu [2], and Kim [3]. Lee [4] Gallaspy [5], and Suprapto [6]; 
Thanh [7], (b) controlling balance by centrifugal force as in the study of Tanaka and Murakami [8], and (c) 
controlling balance by changing the center of gravity as Lee and Ham's research [9]. Among these three 
methods, control of balance using the flywheel has the advantage of being responsive and can be balanced even 
when the vehicle is not moving. 

In two-wheeled robot models that control the balance by using the flywheel, two-wheeled bicycle uses 
the flywheel according to the principle of gyroscope [1, 5-7] to create a balanced torque for the wheels. The 
momentum usually revolves at high speed, so the flywheel dissipates a large amount of energy. This problem 
leads to a difficulty in the operation power of the robot as the robot is only powered by a battery with a limited 
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capacity. In contrast, the two-wheeled bicycle model uses the flywheel according to the principle of inverted 
pendulum [2-4], to create a balanced torque for the car, the flywheel rotates only at very small speed, so the 
energy dissipated by the flywheel is low. Due to the reason, this model is suitable in terms of energy saving 
for the car. Therefore, the authors proposed the self-balancing two-wheeled robot using the flywheel based on 
the principle of inverted pendulum. 

Because two-wheeled bicycles often have to work in different conditions, the carrying capacity may 
vary, the external forces acting on the vehicles may change. It is difficult to find the model of self-balancing 
two-wheeled bicycle, and Two-wheel bicycle can be considered as indeterminate objects [5]. Several control 
algorithms of two-wheeled bicycle have been proposed such as: nonlinear control by Beznol [1], Lee va 
Ham [9], the compensated design using the orbital approach by Gallaspy [5], PD controller by Surpato [6]. 
Due to the uncertainty of two-wheel model, the robust control method in [7] is the most suitable. However, in 
the robust control design method RH. first introduced by McFarlane and Glover in 1992 [10], the controllers 
usually have a high order (controller level is defined as the denominator). The high order controller introduces 
the disadvantage when we use it to control the bicycle. The program is complex. The calculation time is long, 
so the response of the system is slow. Therefore, reducing the order of the controller while ensuring the quality 
of the controller has a significant meaning in practical applications. In order to reduce the controller order, 
there are 2 methods can be followed: 

The first method: this method selects a fixed structure of the order reduction controller and then 
applies optimal algorithms to find the parameters of the order reduction controller so that the standards of the 
robust control are met. The second method: designing a robust controller for an uncertain object will obtain a 
high-order controller, then perform a high-order controller reduction according to the order reduction 
algorithms to obtain a reduced order controller. 

According to the authors, in the first method, the controller can be a low order controller [7], but two 
optimization problems need to be sovle simultaneously (problems in fiding parameters of the controller and 
robust control). This issue leads to difficuty of this method. The parameter of the low order controller may not 
be found if the chosen controller is not suitable. In the second method, the order reduction problem is an 
independent problem, so it always gives the order reduction result as in [11]. Due to that reason, the second 
method has the advantage over the first method because the low order controller can be found in any senario. 

In this paper, the authors proposed the control method of two-wheeled bicycle using model reduction 
algorithm in two steps as follow: (a) design the RH» controller to control the balance of two-wheeled bicycle, 
the found controller is called a full-level controller, and (b) applying order reduction algorithm to reduce order 
of RH» controller to lower order controller while ensuring quality. This step reduction is meant to reduce the 
system response time. 


2. DYNAMIC MODEL AND MATHEMATICAL MODEL OF THE SELF-BALANCING TWO- 
WHEELED BICYCLE 
2.1. Dynamic model of the self-balancing two-wheeled bicycle 

The two-wheeled bicycle model is developed based on the principle of balance using flywheel 
according to the principle of inverted pendulum [2-4]. It is briefly described the principle of balancing of the 
vehicle as follows: if no external torque (torque) is applied to an object or system (or the total torque applied 
to an object is zero), then the total torque of the object will be preserved. 

The vehicle moving by 2 wheels, when the vehicle deviates from the balance position (corresponding 
to a q angle according to vertical axis). The gravity of the vehicle creates a torque that makes the car tend to 
fall down. To maintain a state of equilibrium, we put on the vehicle a flywheel that operates on the principle 
of "the inverted pendulum". This flywheel will rotate around the axis (with an angular acceleration of a) and 
create a torque to compensate the torque generated by the vehicle's gravity. To control the acceleration of the 
flywheel, we uses a DC dc motor with the voltage applied to the motor being U. Then, the problem of balancing 
control becomes the problem of controlling the @ angle (output) by controlling the voltage U (input) applying 
to the motor. The problem requires that the @ angle (output) always go to zero. The self-balancing two-wheeled 
bicycle that the authors built is shown in Figure 1. 

The model machanical parameters: long: 1.19 m; height: 0.5 m; width: 0.4 m; the flywheel weight: 
3.976 kg, diameter: 0.26m; Driving the flywheel using DC motor: 100W-15V-3400 rpm with H-bridge 
driver; Measuring the flywheel velocity by Encoder Sharo 100 pulse; Measuring the q amgle by sensor 
GY-521 MPU-6050; Forward and reverse system consists of a DC motor, H-bridge driver and a remote 
controller. The hardware system is connected to Ardruino microprocessor according to the following block 
diagram as shown in Figure 2. 
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Angle sensor 
GY-521 
MPU-6050 


Figure 2. Schematic structure of bicycle controller 


2.1. Mathematical model of the self-balancing two-wheeled bicycle 

Dynamic model of the self-balancing two-wheeled bicycle is shown in Figure 3. Where: m is the 
bicycle weight (including DC motor), mz is the flywheel weight, h; is the height of the center gravity of the 
bicycle (excluding the flywheel), hz is the height of the center gravity of the flywheel, I; is the inertia torque of 
the bicycle, Iz is the inertia torque of the flywheel, q is the tilt angle of the bicycle corresponding to the vertical 
axis, j is the rotation angle of the flywheel. We have: the absolute velocity of point A is |v,| = h,6. 
The absolute velocity of point B is |vg| = h20. In [5], the author used Lagrange equation to develop the 
dynamic model of the vehicle. 
d (ôT aT ðv 
a T (1) 
where: T is the total kinetic energy of the system, V is the total potential energy of the system, Q; is the external 
force, qi is the generalized coordinate. 











Figure 3. Self-balancing two-wheel bicycle model 


The total kinetic energy of the system defined by: T = T, + T2. Tı, which is the kinetic energy of the 
two-wheeled vehicles, is determined by the following formula: 


1 T A 
T, = 5 malvel* + 587 
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T2, which is the flywheel kinetic energy, is determined by the following formula: 

T, = mlvl? +=1(6 + @) 
We have: 

T = m, val? +5 mylvgl? +56? +216? + hg? + hod (2) 

>T =< (mihi + mA +h + 2)6? +h? + hod (3) 
The total potential energy of the system: 

V = g.cos@.(mh, + mzhz) (4) 
With qi = q, taking (1-4), we get: 

(m,h? + mh? +h +1,)6 + Lö -— g.sin@.(mh, + mhz) = 0 (5) 
With qi =j, taking (1-4), we get: 

LÖ + LÖ = Tp. (6) 
With Tm is the motor shaft torqe. 


Considering a DC de motor with a gear ratio of a:1, the torque of the DC motor driving the flywheel 
is as follows: 


Tm = AKmi = aK, |<", (7) 


with Km is the motor torque constant, Ke is the back-emf constant, R is the resistance of the motor. Substitute 
(7) into (6), we get: 


A 5 U-Keg 
BË + hË = Ty = aKm |$] (8) 
In (5) and (8) are the dynamic system equation. It is clear that the system is nonlinear. Linearizing the model 
and turn it into a state space model. Assume that when the vehicle is operating, the vehicle's inclination angle 


is very small (9 < 10°). Linearizing in (5) around the equilibrium point (9 = g = 0, sin 0 = 0), we have: 


(m,h? + m,h3 + L + L)Ö + Lö — g.9.(m,h, + m,h,) = 0 (9) 





a 5 U-Ked 
hö + hË = Tm = aKm | k 2| (10) 


TakingA, = (m,h? + m,h3 + iL + L); Bı = (mhi + mhz) 


0 =x 
Taking x = | 6 = x, |, is state variable, y = 0,u = U 


P = X3 


We have the state space model describing the system as follow: 








x = Ax + Bu (11) 
y = Cx + Du 
with: 
0 1 0 0 

Big aKmKe a aKm 

A =| (1-2) R(41-12) ;sB= R(A1-I2) C=[1 0 0];D=[0] 
_~ 79 9g sgk K,—“41 — aK,, —41 — 
(41-12) mM e 12R(A1—I2) ™ hR(A1-12) 
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The nominal parameters of the two-wheeled bicycle model are shown in Table 1 as follows: 
explaining research chronological, including research design, research procedure (in the form of algorithms, 
Pseudocode or other), how to test and data acquisition [1-3]. The description of the course of research should 
be supported references, so the explanation can be accepted scientifically [2, 4]. Tables and figures are 
presented center, as shown below and cited in the manuscript. Substituting for the system of (11), we obtain 
the following parameters: 











0 1 0 0 
A =| 47.2048 0 0.0100 |; B =]-0.2230|;C =[1 0 Ol. 
-47.2048 0 -0.1248 2.8541 
Convert the vehicle model into transfer function: 
os) _ —0.223s 
S(s) = U(s)  s3+0.1284s2—47.2s-5.589 (12) 


Table 1. The parameters of the two-wheeled bicycle mode 








Parameter _ Value Unit 
A 0.1105 Kg.m? 
hy 0.105 m 
L 0.03289 Kg.m? 
hy 0.205 m 

mı 10.024 Kg 
mz 3.976 Kg 
Ko 0.045 V.s/Rad 
S 0.045 Nm/A 
R 0.52 Q 
a 1:1 
g 9.81 m/s? 





Remark on two-wheel drive models. The self-balancing two-wheeled bicycle model shows that some 
parameters of self-balancing two-wheeled bicycle are uncertain such as: the changing load volume (leading to 
a change in the center of gravity of the car), the inertia torque of the bicycle changed. Additionally, operating 
two-wheeled bicycle may be influenced by external uncertainties such as: the external force and uncertain noise 
due to the changing of topography. Therefore, The two-wheeled bicycle is the uncertain object. In particular, 
the authors pay the most attention to the uncertainty due to the change of load weight. Specifically, the authors 
consider 4 cases of two-wheeled bicycle carrying different loads as shown in the Table 2. 

Uncertain factors may reduce the accuracy of two-wheeled mathematical models. Therefore, the 
control quality is reduced and the system can even become unstable. Due to the uncertain properties, the various 
control algorithm for the two wheeled bicycle has been proposed: nonlinear control by Beznol [1], Lee và 
Ham [4], the compensated design using the orbital approach by Gallaspy [5], PD controller by Surpato [8]. The 
most suitable algorithm to control the uncertain object was the algorithm in [10]. 


Table 2. Parameters of the two-wheeled bicycle model as the load is different 
Case Load volume Me (kg) Height of the center of gravity hi (m) Moment of inertial, (Kg.m’) 








1 5 0.205 0.6314 
2 5 0.155 0.3609 
3 7 0.055 0.0515 
4 7 0.155 0.409 





3. OPTIMAL DESIGN RH» FOR BALANCE WHEEL PROBLEM 

The structure of the balancing control system for self-balancing two-wheeled bicycle is shown in 
Figure 4. The balancing control system consists of 3 loop controls , namely, loop control the rotation angle of 
the flywheel, loop control the velocity tilt angle of bicycle and loop control the tilt angle of bicylce. The robust 
controller R(s) is used in loop control the tilt angle of bicylce. To design a robust control system for 
self-balancing two-wheeled bicycle, the control structure diagram shown in Figure 4 is used by the authors. 
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w F 7 Two-wheel 
bicycle model 
S(s) 





Figure 4. The robust control structure for self-balancing two-wheeled bicycle 


3.1. Developing the self-balancing two-wheeled bicycle model S, (S) 


Assuming that when the vehicle is in operation, the inclination of the bicycle is very small, we 


linearize in (5) around the equilibrium point (@ = g = 0,sin 0 = 0). We have: 
(m,h?2 + mzh3 + 1, + L)Ö + hö — g.0.(m,h, + mhz) = 0 


LÖ + LÖ = T, = aK, [Severe 


R 


Taking: A, = (mh? + m,h3 + iL + L); Bı = (mh; F mhz) 


0 =x 
Takingx = | = x; |, is state variable, y = 0,u = U* 
$ = x3 


We have the state space model describing the system as follow: 
x = Ax + Bu 
y = Cx + Du 


The system parameters: 











0 1 0 
Big aKmK2 aKm(KetK1) 
A=] (41-1) R(A1-12) R(A1-12) ;B = 
Big Ai At 
(41-12) Knk I2R(A1-12) aRm(Ke + Ka) I2R(A1-12) 


C=[1 0 o0]; D= [0]. 





(13) 
(14) 
(15) 
0 
__aKm 
R(Ay-I2) |; 
A1 


ak,, ————- 
™ 1) R(A1—I2) 


Chosing K, = 2,Kz = 5. Substituting the parameters in Table 1 into (15), the model is converted to the transfer 


function form: 


ecs) __ -0.223s 
U(s)  $3+4.722s2-47.2s—254 


Sm(S) = 





(16) 


To design a robust controller for self-balancing two-wheeled bicycle, the authors followed the steps of 


designing a robust controller RH. according to [10, 12]. We get the robust controller: 


R(s) = 42 


D(s) 


with 


(17) 


H(s) = —2.23.1077s°° — 4.67.10-*5?° — 0.266s?8 — 22.96s?” — 1006s?° — 2.853.10452° 


—5.837.10°s?* — 4.199,1011518 — 9.144,10°s?? — 1.139.108s?? — 1.158.10°s?1 — 9.776.10°s?° 
—6.949,.101%s19 — 2.172.1012s1” — 9.663.101?s16 — 3.71.1013s15 — 1.231.1014s14 
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—3.53.1014513 — 8.74.101451% — 1.862.1015s1* — 3.398,1015s?° — 5.276.101°s? — 6.903.101°s® 
—7.511.10'5s7 — 6.676.101°s°® — 4.721.1075s° — 2.556.1015s* — 9.953.1014s3 — 2.482.1014s? 
—2.977.1013s — 0.00439 

D(s) = 4.971.107 14539 + 2.032.107 19s”? + 2.663.1077s578 + 1.221.107 *5?” + 9.72.107 357° 
+0.3918s?5 + 10.14574 + 187.1573 + 261257? + 2.862.104s7! + 1.088.107s51® + 2.523.10°s?° 
+1.82.10°s?? + 5.428.107s17 + 2.273.108s16 + 8.005.108s1® + 2.372.10°s1* + 5.9.10°s?3 
+1.225.1029s?2 + 2.107.101°s14 + 2.962.1029s?° + 3.341.101°s? + 2.941.1019s® 
+1.931.107°s” + 8.743.10°s® + 2.286.109s° + 1.519.108s* — 5.226.107s? + 3.6.10 °s? 


+5.32.107??s 


3.2. Compare the robust controller with another controller 

The bicycle with the varying parameter is controlled by PID control method. The result is then used 
to compare to the case which the robust controller is applied. Simulation diagram of self-balancing 
two-wheeled bicycle control system using robust controller and PID controller are shown in Figure 5. 
Simulation results of self-balancing two-wheeled bicycle control system when the parameters of model are 
rated and when the model parameters change, Initially, the bicycle deviates 60 = 5 (rad) from the vertical 


axis. Parameters of PID controller: Kp = -450, Kı = -30, Kp = -15. The results shown in Figure 6. 






num(s) 
den(s) 
30th-order robust controller. Saturation’ 












Two-wheel bicycle model. 








Step. 











PID Controller. Saturation 1 


Two-wheel bicycle model., 


Figure 5. Simulation diagram of self-balancing two-wheeled bicycle control system using robust controller 
and PID controller 


Remark by the simulation result in both cases the nominal parameters of the bicycle and the variable 
parameters of the bicycle due to the load varied. PID controller can only balance the bicycle as the bicycle 
parameters is niminal and in case 3. PID controller has not worked in the case 1, 2, and 4. The robust controller 
did work in all 4 cases. It can be seen that the robust controller was able to balance the system even the system 
parameters are varied (load and the height of the center of gravity of the bicycle). The robust controller has the 
advantage over the PID controller 

Howerver, the 30" order controller could lead to the difficulty of operating the balancing process. Due 
to the complex program, the long processing time, the low system response, the system will not be able to adapt 
the requirements of real-time applications and can become unstable. For that reason, reducing the order of the 
controller is needed to simplify the program. The system response is thereby increased, while the robustness is 
ensured. 
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Figure 6. The system output response of self-balancing two-wheeled bicycle control system using robust 
controller and PID controller 


4. STOCHASTIC BALANCE TRUNCATION ALGORITHM BASED ON SCHUR ANALYSIS 


4.1. Model reduction problem 
Given a linear, continuous, time-invariant, MIMO system described by the following state space 


model: 
x = Ax + Bu y = Cx (18) 


where, x € R”,u E€ R?,y € R1, A € R™", B e R™P,C e R™". The goal of the order reduction problem for 
the model described by state space model given in (17) is to find the model described by state space model: 


X, = A,X, + Bu 
Yr = CX, (19) 


where, x, E R’,u, E RP, y, E R1, A4, E R", B, E R™,C, E RU voir «Kn. So that the model described by 
in (19) can replace the model described by the (18) applications in analysis, design, and control of the system. 
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4.2. Stochastic balanced truncation algorithm based on Schur analysis 

Most of the model reduction algorithms have published in the world only apply to stable high order 
linear models (the roots of the charateristic equation are negative) [13-15]. However, many high order 
mathematical models are unstable in reality such as the model in section 3. Therefore, the order reduction 
algorithm should be applicable to reduce the order for the unstable linear system. There are two basic methods 
for model reduction of unstable system. The first method (indirect order reduction algorithm). This algorithm 
divides the unstable original system into stable and unstable components, then applies the order reduction 
algorithm to the stable components [16-24]. At the end, to get the order of reduction of the root system, we add 
the reduced stable components with the unstable components. 

The second method (direct order reduction algorithm). This algorithm modifies and adjusts the order 
reduction algorithms so that these algorithms can perform order reduction regardless of whether the original 
system is stable or unstable [25-29]. In the content of this paper, the author introduces the stochastic balanced 
truncation algorithm based on Schur analysis [23, 24]. This is a order reduction algorithm applied to the 
unstable system by indirect order reduction method. The specific contents of the algorithm are as follows: 
Input: The system (A, B, C) (stable or unstable) described in (18) has a representation of the form of the transfer 
function: G(s): = C(sI — A)71B. 

Step 1: Find the controllability grammian P and observability grammian Q by solving the following Lyapunov 
and Riccati equations: 


AP + PAT + BB" = 0; By = PC’ + BD"; QA + ATQ + (QBy —C™)(—DD")(QBy — C7)’ =0 


Step 2: Find the Schur decomposition for PQ in both ascending and descending order, respectively, 


Ags age Hee 

0 0 A 

Step 3: Find the left/right orthonormal eigen-bases of PQ associated with the k" big Hankel singular values of 
the all-pass phase matrix (W * (s))~*G(s). 


Ay one one 
Op ace | VEPQV, = 


VI PQV, = 7 
0 0 An 








k 


Vp = 


ge 
V, = VR,SMALL» VL BIG 











k 
oo 
Va Bic , VR SMALL 


Step 4: Find the SVD of (VigicVrpic) = ULV 
Step 5: Form the left/right transformation for the final k® order reduced model. 


SL BIG = V gıGUZ (1: k, 1: kl; SR BIG = Vr pigV2(1: k, 1: kj)? 


Step 6: Calculate (A,, B,,C,) = (SP argASpia SLBrcB, CSR p16 )- 
Output: The reduced system (A,,, B,, C,). 


5. APPLIED LQG ALGORITHM FOR ROBUST CONTROL PROBLEM OF TWO-WHEELED 
BICYCLE 
5.1. The reduced controller of sef-balancing two-wheeled bicycle 
The full order RH.. controller is designed as (17), which is a 30th order controller. To obtain low 
controller, we perform order reduction of RH.. controller in accordance with the stochastic balanced truncation 
algorithm based on Schur analysis in section 4. The results of the order reduction controller are shown in 
Table 3. 


5.2. Controlling the two-wheeled bicycle using the reduced 4th and 5th order controller 

Using the reduced 5* order controller in Table 3 controls the balancing system for two-wheeled 
bicycle having the model as (16). The performance is compared to the performance of the original (30" order) 
controller. The simulation diagram of two-wheeled bicycle system using the original controller and reduced 
controllers is shown in Figure 7. 
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Table 3. Results of model reduction of the 30" order robust controller 
Order Transfer function — R,(s) 
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s>+2006s4+1.135.10453-1193s7+6.179.107 135—4.24.10- 14 
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s4+2000s3—206.552+2.369.10~14s—3.026.107-15 
—4.485.10°s?+4.047.105s*—9.163.104s—5.881.10* 


s3—0,1032s7+5.41,10717s—3,669,10718 








4 





3 









0th-order robust controller. 


-4,485e0859-5 241207s*-3.026e08s°-8, 389e08s2-1.17 7209-6. 796208 
s5+2008s*+1.138e4s9-1193s2+6. 1792-1954. 240-14 


Ath-orde robust cont oller (Bstmr) Setwaton24 


Twowheel bicycle model 





Scope 


-4,485e06s4-2 655207s2-1.141208s2-1.833e08s-1.176208 
s4+2000s?.206.8s2+1.082e-13s-7.338e-15 






4thorder robust controller (Bstmr) 


Two-wheel bicycle model. 





O -4.485208s9+4,047208s2.9.163e45-5.88164 
$7-0,1032s2+5,416-175-3,669e-18 


Zd- order robust controller (Bstmr) Saturation31 





Two-wheel bicycle model... 





Figure 7. The simulation diagram of two-wheeled bicycle system using the original controller and reduced 
controllers 


The simulation results of self-balancing two-wheeled bicycle control system when the parameters 


model are nominal and the bicycle is initially deviated 0 = 1/180 (rad) from the vertical axis. The results 
shown in Figure 8. The simulation results of self-balancing two-wheeled bicycle control system when 
the vehicle is loaded and the bicycle is initially deviated 0 = 77/180 (rad) from the vertical axis. The results 


shown in Figure 9. 
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Figure 8. Output response of the self-balancing two-wheeled bicycle control system in no load condition 
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Figure 9. Output response of the self-balancing two-wheeled bicycle control system when the system is 


R,(s) = 


loaded condition (continue) 


Compared the result between the system using the original controller and the system using the reduced 
controller according to the balanced truncation and other algorithms. The author compared the stochastic 
balanced truncation algorithm based on Schur analysis with the balanced truncation algorithm proposed by 


Moore [30]. This is the most commonly used order reduction algorithm. In Matlab, the command balancmr is 
used to perform the balanced truncation. We get the reduced 4" order controller. 


—4.485.10°s*—5.351.108s3+7.513.107s7+2.822.107s+1.307.107 


s*+2000s3—206.5s24+1.258.10-19s—4.767.10~ 12 
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The simulation is performed with the nominal parameters of the two-wheeled bicycle model and initially 
deviated from the vertical by an angle 0 = oT (rad). The result shown in Figure 10. 
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Figure 10. Output response of the self-balancing two-wheeled bicycle control system using the reduced 
4" order controller 


5.3. Evaluated results 

The reduced 4" and 5" order controller arcoording to the stochastic balanced truncation algorithm 
based on Schur analysis can be used to control the two-wheeled bicycle model. The output response of the 
reduced 5" order controller is almost identical to the output response of original controller. The output response 
of the reduced 4" order controller is different from that of the original controller. The reduced 3" order 
controller cannot control the balance of the two-wheeled bicycle model. 

Compared the result of the two-wheel balancing control system between the system using the reduced 
controller according to the balanced truncation algorithm based on Schur analysis and the system using the 
balanced truncation algorithm (balancmr): We see that the control system using the reduced 4th order controller 
according to the balanced truncation algorithm based on Schur analysis ensure the stable balance of the two- 
wheeled bicycle when the bicycle deviates from vertical and when the parameters of the model change, while 
the control system using the reduced 4th order controller according to the balanced truncation algorithm 
(balancmr) does not. 


6. CONCLUSION 

The paper has developed, modeled a two-wheeled self-balancing bicycle model and designed a robust 
controller to control the balance of two-wheeled bicyle. The paper also introduces the stochastic balanced 
truncation algorithm based on Schur analysis and applies this algorithm to reduce the high order robust 
controller using to control the balance of two-wheeled bicyle. In particularly, the reduced 4" and 5" order 
controller can replace the original controller (30"-order) while the performance of the control system is 
ensured. Using the reduced controller simplify the program, so the computational time is reduced. Therefore, 
the system respose is improved, and the requirements in real-time application are met. The simulation results 
show the correctness of the model reduction algorithm and the robust control algorithm of two-wheeled 
self-balancing two-wheeled bicycle. 
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